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Abstract. We develop a tropical analogue of the classical double description method 
allowing one to compute an internal representation (in terms of vertices) of a polyhedron 
defined externally (by inequalities). The heart of the tropical algorithm is a characteri- 
zation of the extreme points of a polyhedron in terms of a system of constraints which 
define it. We show that checking the extremality of a point reduces to checking whether 
there is only one minimal strongly connected component in an hypergraph. The latter 
problem can be solved in almost linear time, which allows us to eliminate quickly redun- 
dant generators. We report extensive tests (including benchmarks from an application to 
static analysis) showing that the method outperforms experimentally the previous ones by 
orders of magnitude. The present tools also lead to worst case bounds which improve the 
ones provided by previous methods. 



Introduction 

Tropical polyhedra are the analogues of convex polyhedra in tropical algebra. The latter 
deals with structures like the max-plus semiring ]R max (also called max-plus algcbvci^ which 
is the set R U {— oo}, equipped with the addition x © y := max(x, y) and the multiplication 
x ® y := x + y. 

The study of the analogues of convex sets in tropical or max-plus algebra is an active 
research topic, and has been treated under various guises. It arose in the work of Zim- 
merman [Zim77j, following a way opened by Vorobyev [Vor67] , motivated by optimization 
theory. Max-plus cones were studied by Cuninghame-Green |CG79| . Their theory was 
independently developed by Litvinov, Maslov and Shpiz [LMSOlJ (see also [MS92| ) with 
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motivations from variations calculus and asymptotic analysis, and by Cohen, Gaubert, and 
Quadrat |CGQ04| who initiated a "geometric approach" of discrete event systems |CG Q99|, 
further developed in [Kat07, DLGKL09J. Other motivations arise from abstract convexity, 
see the book by Singer jSin97| . and also the work of Briec and Horvath [BH04j . The field 
has attracted recently more attention after the work of Develin and Sturmfels [DS04 ] , who 
pointed out connections with tropical geometry, leading to several works by Joswig, Yu, 
and the same authors |Jos05l IDY071 IJSY071 lJosT)9] . 

A tropical polyhedron can be represented in two different ways, either internally, in 
terms of extreme points and rays, or externally, in terms of linear inequalities (see Sect. Q] 
for details). As in the classical case, passing from the external description of a polyhedron 
to its internal description is a fundamental computational issue. This is the object of the 
present paper. 

Butkovic and Hegedus [BH84J gave an algorithm to compute the generators of a tropical 
polyhedral cone described by linear inequalities. Gaubert gave a similar one and derived the 
equivalence between the internal and external representations [Gau92, Ch. Ill] (see [GK09 
for a recent discussion). Both algorithms rely on a successive elimination of inequalities, 
but have the inconvenience of squaring at each step the number of candidate generators, 
unless an elimination technique is used, as in the Maxplus toolbox of Scilab |CGMQ . 
Joswig developed a different approach, implemented in Polymake [GJj . in which a tropical 
polytope is represented as a polyhedral complex [DS04, Jos09]. 

The present work grew out from two applications: to discrete event systems |Kat07t 
DLGKL09J, and to software verification by static analysis |AGG08| . In these applications, 
passing from the external to the internal representation is a central difficulty. A further 
motivation originates from mean payoff games [AGG09b]. These motivations are reviewed 
in Section [21 

Contributions. We develop a new algorithm which computes the extreme elements of trop- 
ical polyhedra. It is based on a successive elimination of inequalities, and a result (Th. 14. ip 
allowing one, given a polyhedron V and a tropical halfspace H, to construct a list of can- 
didates for the generators of V H H. The key ingredient is a combinatorial characterization 
of the extreme generators of a polyhedron defined externally (Th. 13.5 1 and 13.7] ): we reduce 
the verification of the extremality of a candidate to the existence of a strongly connected 
component reachable from any other in a directed hypergraph. We include a complexity 
analysis and experimental results (Sect.HJ, showing that the new algorithm outperforms the 
earlier ones, allowing us to solve instances which were previously by far inaccessible. Our 
result also leads to worst case bounds improving the ones of previously known algorithms. 

1. Definitions: tropical polyhedra and polyhedral cones 

The neutral elements for the addition and multiplication (g), i.e., the zero and the 
unit, will be denoted by := — oo and 1 := 0, respectively. The tropical analogues of the 
operations on vectors and matrices are defined naturally. The elements of M^axi the dth 
fold Cartesian product of M max , will be thought of as vectors, and denoted by bold symbols, 
like x = (jci, . . . ,x d ). 

A tropical halfspace is a set of the vectors x = (xi) € H^ax verifying an inequality 
constraint of the form 



max a, + x ; < max bi + Xi, ai,bi G M, 

l<i<d l<i<d 
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A tropical polyhedral cone is defined as the intersection of n halfspaces. It can be equivalently 
written as the set of the solutions of a system of inequality constraints Ax < Bx. Here, 
A = (a,ij) and B = (by) are n x d matrices with entries in M max , concatenation denotes 
the matrix product (with the laws of M m ax)i an d < denotes the standard partial ordering 
of vectors. For sake of readability, tropical polyhedral cones will be simply referred to as 
polyhedral cones or cones. 

Tropical polyhedral cones are known to be generated by their extreme rays [GK06, 
IGK071 IBSS07] . Recall that a ray is the set of scalar multiples of a non-zero vector u. It 
is extreme in a cone C if u G C and if it = v © w with v,w £ C implies that u = v or 
u = w. A finite set G = (g r )i£i of vectors is said to generate a polyhedral cone C if each 
g l belongs to C, and if every vector x of C can be written as a tropical linear combination 
®i^i9 % of the vectors of G (with Aj 6 M max ). Note that in tropical linear combinations, 
the requirement that Aj be nonnegative is omitted. Indeed, = — oo < A holds for all scalar 

A G Imax- 

The tropical analogue of the Minkowski theorem GK07, BS S07] shows in particular 
that every generating set of a cone that is minimal for inclusion is obtained by selecting 
precisely one (non-zero) element in each extreme ray. 

A tropical polyhedron of M max is the affine analogue of a tropical polyhedral cone. It is 
defined by a system of inequalities of the form Ax © e < Bx © d. It can be also expressed 
as the set of the tropical affine combinations of its generators. The latter are of the form 
©j 6 2 Aji>* © © je j fJ-jf J , where the (v l )i<zi are the extreme points, the (r J )j<zj a set formed 
by one element of each extreme ray, and © i Aj = 1. It is known |CGQ04, GK07] that 



every tropical polyhedron of M max can be represented by a tropical polyhedral cone of 
thanks to an analogue of the homogenization method used in the classical case (see [Zie98, 
Sect. 1.5]). Then, the extreme rays of the cone are in one-to-one correspondence with the 
extreme generators of the polyhedron. That is why, in the present paper, we will only state 
the main results for cones, leaving to the reader the derivation of the affine analogues, along 
the lines of [GK07| . 

In the sequel, we will illustrate our results on the polyhedral cone C given in Fig. [H 
defined by the system in the right side. The left side is a representation of C in barycen- 
tric coordinates: each element (xi,X2,xs) is represented as a barycenter with weights 
(e Xl , e X2 , e X3 ) of the three vertices of the outermost triangle. Then two elements of a same 
ray are represented by the same point. The cone C is depicted in solid gray (the black bor- 
der is included), and is generated by the extreme elements g° = (0,0,0), g 1 = (—2,1,0), 
<? 2 = (2,2,0), and g 3 = (0,0,0). 

2. Motivations from static analysis, discrete event systems, and mean pay- 
off games 

Tropical polyhedra have been recently involved in static analysis by abstract interpreta- 
tion [AGG08J. It has been shown that they allow to automatically compute complex invari- 
ants involving the operators min and max which hold over the variables of a program. Such 
invariants are disjunctive, while most existing techniques in abstract interpretation are only 
able to express conjunctions of affine constraints, see in particular [CC77[ ICH781 iMinOl] . 

For instance, tropical polyhedra can handle notorious problems in verification of mem- 
ory manipulations. Consider the well-known memory string manipulating function memcpy 
in C. A call to memcpy(dst, src,n) copies exactly the first n characters of the string buffer 
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Figure 2: memcpy invariant 



src to dst. In program verification, precise invariants over the length of the strings are 
needed to ensure the absence of string buffer overflows. Recall that the length of a string 
is defined as the position of the first null character in the string. To precisely analyze the 
function memcpy, two cases have to be distinguished: 

(i) either n is strictly smaller than the source length len_src, so that only non-null char- 
acters are copied into dst, hence len_dst > n, 

(ii) or n > len_src and the null terminal character of src will be copied into dst, thus 
len_dst = len_src. 

Thanks to tropical polyhedra, the invariant min(len_src, n) = min(len_dst, n), or equiva- 
lently max(— len_src, — n) = max(— len_dst, — n), can be automatically inferred. It is the 
exact encoding of the disjunction of the cases (i) and (ii). The invariant is represented by 
the non-convex set of M 3 depicted in Figure [2j In the application to static analysis, the 
performance of the algorithm computing the extreme elements of tropical polyhedra plays 
a crucial role in the scalability of the analyzer (see [AGG08] for further details). 

A second motivation arises from the "geometric approach" of max-plus linear discrete 
event systems |CGQ99|, in which the computation of feedbacks ensuring that the state of 
the system meets a prescribed constraint (for instance that certain waiting times remain 
bounded) reduces [Kat07j to computing the greatest fixed point of an order preserving map 
on the set of tropical polyhedra. Similar computations arise when solving dual observability 
problems [DLGKL09]. Again, the effective handling of these polyhedra turns out to be the 
bottleneck. 

A third motivation arises from the study of mean payoff combinatorial games. In 
particular, it is shown in |AGG09b] that checking whether a given initial state of a mean 
payoff game is winning is equivalent to finding a vector in an associated tropical polyhedral 
cone (with a prescribed finite coordinate). This polyhedron consists of the super- fixed points 
of the dynamic programming operator (potentials), which certify that the game is winning. 
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3. Characterizing extremality from inequality constraints 

3.1. Preliminaries on extremality 

The following lemma, which is a variation on the proof of Th. 3.1 of [GK07] and on 
Th. 14 of [BSS07], shows that extremality can be expressed as a minimality property: 

Proposition 3.1. Given a polyhedral cone C C K max , g is extreme if and only if there 
exists 1 < t < d such that g is a minimal element of the set { x £ C \ x% = g t }, i.e. g € C 
and for each x £ C , x < g and Xt = g t implies x = g. In that case, g is said to be extreme 
of type t. 

In Fig. [3j the light gray area represents the set of the elements (x\, X2,x%) of K max such 
that (x\, X2, X3) < g 2 implies x\ < g\. It clearly contains the whole cone except g 2 , which 
shows that g 2 is extreme of type 1. 

A tropical segment is the set of the tropical linear combinations of two points. Using the 
fact that a tropical segment joining two points of a polyhedral cone C yields a continuous 
path included in C, one can check that g is extreme of type t in C if and only if there is a 
neighborhood N of g such that g is minimal in {x GCniV [ Xt = g t }. Thus, extremality 
is a local property. 

Finally, the extremality of an element g in a cone C can be equivalently established by 
considering the vector formed by its non-0 coordinates. Formally, let supp(a;) := {i \ Xi 7^ 0} 
for any x € K max . Then g is extreme in C if and only if it is extreme in {x € C | supp(x) C 
supp(^)}. This allows to assume that supp(g) = {1, . . . , d} without loss of generality. 

3.2. Expressing extremality using the tangent cone 

For now, the polyhedral cone C is supposed to be defined by a system Ax < Bx of n 
inequalities. 

Consider an element g of the cone C, which we assume, from the previous discussion, 
to satisfy supp(g) = { 1, . . . ,d}. In this context, the tangent cone of C at g is defined as 
the tropical polyhedral cone T(g, C) of K max given by the system of inequalities 

max Xi < max Xj for all k such that A^g = B^g, (3-1) 

ieargmax(A fe g) jearg max(B fc g) 

where for each row vector c G I^max, argmax(cg) is defined as the argument of the maximum 
eg = maxi<j< ( i(cj+gi i ), and where A\~ and B^ denote the kth rows of A and B, respectively. 
The tangent cone T(g,C) provides a local description of the cone C around g: 

Proposition 3.2. There exists a neighborhood N of g such that for all x € N, x belongs 
to C if and only if it is an element of g + T(g, C). 

As an illustration, Fig. H] depicts the set g 2 + T(g 2 ,C) (in semi-transparent light gray) 
when C is the cone given in Fig. [TJ Both clearly coincide in the neighborhood of g 2 . 

Since extremality is a local property, it can be equivalently characterized in terms of 
the tangent cone. Let 1 be the element of K max whose all coordinates are equal to 1. 

Proposition 3.3. The element g is extreme in C iff the vector 1 is extreme in T(g,C). 



52 



XAVIER ALLAMIGEON, STEPHANE GAUBERT, AND ERIC GOUBAULT 



X 3 



Xi X 2 

Figure 3: Extremality of g 2 
x 3 m (0,0,1) 




(0,1,1) 



Figure 5: The {0,1 }-elements of T(g 2 ,C) 




Figure 4: The set g 2 + T(g 2 ,C) 




Figure 6: A directed hypergraph 



The problem is now reduced to the characterization of the extremality of the vector 1 
in a { 0, 1 }-cone, i.e. a polyhedral cone defined by a system of the form Cx < Dx where 
C,D € { 0, 1 } nxd . The following proposition states that only { 0, 1 }- vectors, i.e. elements 
of the tropical regular cube {0,1 } d , have to be considered: 



Proposition 3.4. Let T> C 



^max 



be a { 0, 1 }-cone. Then 1 is extreme of type t if and only 



if it is the unique element x of T> n { 0, 1 } satisfying x% = 1. 

The following criterion of extremality is a direct consequence of Prop. 13.31 and [37 

Theorem 3.5. Let C C K^ax be a polyhedral cone. Then g € C is extreme of type t if and 
only if the vector 1 is the unique {§),!} -element of the tangent cone T(g,C) whose t-th 
coordinate is 1. 

Figure [5] shows that in our running example, the { 0, 1 }-elements of T(g 2 ,C) distinct 
from 1 (in squares) all satisfy X\ = 0. Naturally, testing, by exploration, whether the set 
of 2 d ~ l {0,1 }-elements x verifying Xt = 1 belonging to T(gr, C) consists only of 1 does not 
have an acceptable complexity. Instead, the approach of the next section will rely on the 
equivalent formulation of the criterion of Th. I3.5t 



VZe {l,...,d}, [Vx£T(g,C)n{®,t} d , Xl 



x t = . 



(3.2) 



3.3. Characterizing extremality with directed hypergraphs 

A directed hypergraph is a couple (N, E) such that each element of E is of the form 
(T, H) with T,H CN. 
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The elements of N and E are respectively called nodes and hyperedges. Given a hyper- 
edge e = (T, H) E E, the sets T and H represent the tail and the head of e respectively, 
and are also denoted by T(e) and H(e). Figure [6] depicts an example of hypergraph whose 
nodes are u, v, w, x, y, t, and of hyperedges e\ = ({it}, {v}), e2 = ({v}, {w}), e% = ({w}, {u}), 
e 4 = ({v,w},{x,y}), and e 5 = ({to, y}, {*}). 

Reachability is extended from digraphs to directed hypergraphs by the following recur- 
sive definition: given u, v E N, then v is reachable from u in Ti, which is denoted u ~~^>u v, 
if one of the two conditions holds: u = v, or there exists e E E such that v E H(e) and all 
the elements of T(e) are reachable from u. In our example, t is reachable from ?/. 

The size size(^) of a hypergraph U = (N, E) is defined as \N\ + J2 eeE (\ T ( e )\ + \ H ( e )\)- 
In the rest of the paper, directed hypergraphs will be simply referred to as hypergraphs. 

We associate to the tangent cone T(g,C) the hypergraph 7i(g,C) = (N,E) defined by: 

N = {l,...,d} E = { (argmax( J B fc si),argmax(^ fc gr)) | A k g = B k g, 1 < k < n } . 

The extremality criterion of Eq. (13, 2p suggests to evaluate, given an element of T(g, C) n 
{0,1 } d , the effect of setting its l-th coordinate to the other coordinates. Suppose that it 
has been discovered that x\ = implies Xj t = ■ ■ ■ = Xj n = 0. For any hyperedge e of 
H(g,C) such that T(e) C ■ ■ ■ ,j n }i x satisfies: max ie ^( e j Xi < maXj GT r e \Xj = 0, so 

that Xi = for all % E H{e). Thus, the propagation of the value from the Z-th coordinate 
to other coordinates mimicks the inductive definition of the reachability relation from the 
node I in H(g,C): 

Proposition 3.6. For all I E { 1, . . . ,d}, the statement given between brackets in Eq. (j3.2|) 

holds if and only if t is reachable from I in the hypergraph TL(g,C). 

Hence, the extremality criterion can be restated thanks to some considerations on the 
strongly connected components of 7i(g,C). The strongly connected components (Sees for 
short) of a hypergraph Ti are the equivalence classes of the equivalence relation defined 
by u =u v if u -w^ v and v u. They form a partition of the set of nodes of H. They 
can be partially ordered by the relation defined by C\ <u C2 if C% and C2 admit a 
representative u and v respectively such that v u (beware of the order of v and u in 
v u). Then Prop. [3T61 and Th. 13.51 imply the following statement: 

Theorem 3.7. Let C C M^ ax be a polyhedral cone, and g E C. Then g is extreme if and 
only if the set of the Sees of the hypergraph 7i(g,C), partially ordered by ^^( S) c) ; admits a 
least element. 

This theorem is reminiscent of a classical result, showing that a point of a polyhedron 
defined by inequalities is extreme if and only if the family of gradients of active inequalities 
at this point is of full rank. Here, the hypergraph encodes precisely the subdifferentials 
(set of generalized gradients) of the active inequalities but a major difference is that the 
rank condition must be replaced by the above minimality condition, which is essentially 
stronger. Indeed, using this theorem, it is shown in [AGK09] that an important class of 
tropical polyhedra has fewer extreme rays than its classical analogue. 

An algorithm due to Gallo et al. [GLPN93] shows that one can compute the set of nodes 
that are reachable from a given node in linear time in an hypergraph. The following result 
shows that one can in fact compute the minimal Sees with almost the same complexity. 
The algorithm is included in the extended version of the present paper |AGG09c] . Although 
it shows some analogy with the classical Tarjan algorithm, the hypergraph case differs 
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critically from the graph case in that one cannot compute all the Sees using the same 
technique. 

Theorem 3.8. The set of minimal Sees of a hypergraph T~L = (N, E) can be computed in 
time 0(s\ze(H) x a(|iV|)), where a denotes the inverse of the Ackermann function. 



4. The tropical double description method 

Our algorithm is based on a successive elimination of inequalities. Given a polyhedral 
cone C defined by a system of n constraints, the algorithm computes by induction on k 
(0 < k < n) a generating set Gk of the intermediate cone defined by the first k constraints. 
Then G n forms a generating set of the cone C. Passing from the set Gk to the set Gk+i relies 
on a result which, given a polyhedral cone fC and a tropical halfspace H = { x \ ax < bx }, 
allows to build a generating set G' of K, fl H from a generating set G of /C: 

Theorem 4.1. Let K, be a polyhedral cone generated by a set G C and Tl = { x \ 

ax < bx} a tropical halfspace (a, b € IR^*^- Then the polyhedral cone K. fl H is generated 
by the set {g G G \ ag < bg} U { {ah)g © (bg)h \ g,h € G, ag < bg, and ah > bh }. 

For instance, consider the cone defined in Fig. Q] and the constraint x<i < X3 + 2.5 
(depicted in semi-transparent gray in Fig. [7J. The three generators g , g 2 , and gi 3 satisfy 
the constraint, while g° does not. Their combinations are the elements h 1 ' , h 2 ' , and ft. 3 ' 
respectively. The resulting algorithm is given in Figure [U As in the classical case, this 
inductive approach produces redundant generators, hence, the heart of the algorithm is the 
extremality test in Line 10. We use here the hypergraph characterization (Theorems 13.71 
and ES]) • 

Complexity analysis. The complexity of the elementary step of ComputeExtreme, i.e. the 
computation of the elements provided by Th. 14.11 and the elimination of non-extreme ones 
(Lines [7] to [T3j) , can be precisely characterized to 0(nda(d) |G| 2 ), where G is the generating 
set of the last intermediate cone. By comparison, for classical polyhedra, the same step in 
the refined double description method by Fukuda and Prodon }FP96| takes a time 0(n |G| 3 ). 
Note that |G| can take values much larger that d. 
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1: 
2: 
3: 
4: 
5: 
6: 
7: 
8: 
9: 
10: 
11: 
12: 
13: 



procedure ComputeExtreme^, B, n) t> A,B g R; 
if n = then > Base case 



return the tropical canonical basis (ei)i<i<d 



else l> Inductive case 



split Ax < Bx into Cx < Dx and ace < bx, with C, -D g ln„ 1 ' x ' ! and a, 6 g R, 
G := ComputeExtreme(C, D, n — 1) 

G- :={9'eG ag' < bg' }, G> :={g J 6G > 6g^ }, // := G- 

for all g i g G^ and gi g G> do 



-max 



h := (ag^g* ® (bg*)gi 



if h is extreme in { x \ Ax < Bx } then 



append nh to H , where n is the opposite of the first non-(D coefficient of h 



end 
done 



14: end 

15: return H 

16: end 

Figure 9: Our main algorithm computing the extreme rays of tropical cones 

The overall complexity of the algorithm ComputeExtreme depends on the size of the 
sets returned in the intermediate steps. In classical geometry, the upper bound theorem of 
McMullen [McM70] shows that the maximal number of extreme points of a convex polytope 
in R rf defined by n inequality constraints is equal to 



The polars of the cyclic polytopes (see |Zie98j ) are known to reach this bound. Allamigeon, 
Gaubert, and Katz [AGK09] showed that a similar bound is valid in the tropical setting. 

Theorem 4.2 ([AGK09J). The number of extreme rays of a tropical cone in (MU {— oo}) d 
defined as the intersection of n tropical half-spaces cannot exceed U(n + d,d — 1) = 0{{n + 



The bound is asymptotically tight for a fixed n, as d tends to infinity, being approached 
by a tropical generalization of the (polar of) the cyclic polytope |AGK09j . The bound is 
believed not to be tight for a fixed d, as n tends to infinity. Finding the optimal bound is an 
open problem. By combining Theorem 14.21 Theorem 13.8^ and Theorem 13.71 we readily get 
the following complexity result, showing that the execution time is smaller in the tropical 
case than in the classical case, even with the refinements of [FP96J. 

Proposition 4.3. The hypergraph implementation of the tropical double description method 
returns the set of extreme rays of a polyhedral cone defined by n inequalities in dimension 
d in time 0(n 2 da(d)G max ), where G max is the maximal number of extreme rays of a cone 
defined by a subsystem of inequalities taken from Ax < Bx. In particular, the time can be 
bounded by 0{n 2 da(d){n + d)^ 1 ). 

Alternative approaches. The existing approachs discussed in the introduction have a struc- 
ture which is similar to ComputeExtreme. However, their implementation in the Maxplus 
toolbox of SciLAB |CGMQ| and in our previous work [AGG08] relies on a much less efficient 
elimination of redundant generators. Its principle is the following: an element h is extreme 
in the cone generated by a given set H if and only if h can not be expressed as the tropical 
linear combination of the elements of H which are not proportional to it. This property can 
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Table 1: Benchmarks on a single core of a 3 GHz Intel Xeon with 3 Gb RAM 





d 


n 


# final 


# inter. 


T (s) 


T' (s) 


T/T' 


rndlOO 


12 


15 


32 


59 


0.24 


6.72 


0.035 




rndlOO 


15 


10 


555 


292 


2.87 


321.78 


8.9 • 10" 


3 


rndlOO 


15 


18 


152 


211 


6.26 


899.21 


7.0 • 10" 


3 


rnd30 


17 


10 


1484 


627 


15.2 


4667.9 


3.3 • 10" 


3 


rndlO 


20 


8 


5153 


1273 


49.8 


50941.9 


9.7 • 10" 


4 


rndlO 


25 


5 


3999 


808 


9.9 


12177.0 


8.1 • 10" 


-1 


rndlO 


25 


10 


32699 


6670 


3015.7 








cyclic 


10 


20 


3296 


887 


25.8 


4957.1 


5.2 • 10" 


3 


cyclic 


15 


7 


2640 


740 


8.1 


1672.2 


5.2 • 10" 


3 


cyclic 


17 


8 


4895 


1589 


44.8 


25861.1 


1.7 • 10" 


3 


cyclic 


20 


8 


28028 


5101 


690 


45 days 


1.8 • 10" 


4 


cyclic 


25 


5 


25025 


1983 


62.6 


8 days 


9.1 • 10" 


5 


cyclic 


30 


5 


61880 


3804 


261 








cyclic 


35 


5 


155040 


7695 


1232.6 











# var 


# lines 


T (s) 


T (s) 


T/T 


oddcvcn8 


17 


118 


7.6 


152.1 


0.050 


oddcvcn9 


19 


214 


128.0 


22101.2 


5.8 ■ 10~ 3 


oddcvcnlO 


21 


240 


1049.0 







be checked in 0{d x \H\) time using residuation (see [BSS07] for algorithmic details). In 
the context of our algorithm, the worst case complexity of the redundandy test is therefore 
0(d |G| 2 ), where G is the set of the extreme rays of the last intermediary cone. This is much 
worse that our method in 0(nda(d)) based on directed hypergraphs, since the cardinality 
of the set G may be exponential in d (see Theorem 14. 2p . This is also confirmed by our 
experiments (see below). 

We next sketch a different method relying on arrangement of tropical hyperplanes 
(arrangements of classical hyperplanes yield naive bounds). Indeed, Theorem 13.71 implies 
that every extreme ray belongs to the intersection of d — 1 tropical hyperplanes, obtained 
by saturating d — 1 inequalities among the n + d taken from Ax < Bx and Xi > — oo, for 
i € [d]. The max-plus Cramer theorem (see [AGG09a] and the references therein) implies 
that for generic values of the matrices A,B, every choice of d — 1 saturated inequalities 
yields at most one candidate to be an extreme ray, which can be computed in 0(d 3 ) time. 
This yields a list of 0{{n + d) d ~ l ) candidates, from which the extreme rays can be extracted 
by using the present hypergraph characterization (Theorems 13.71 and I3.8[) . leading to a 
0((nda(d) + d 3 )(n + c?) d_1 ) execution time, which is better than the one of Proposition 14.31 
by a factor n/a(d) when n d. However, the resulting algorithm is of little practical 
use, since the worst case execution time is essentially always achieved, whereas the double 
description method takes advantage of the fact that G max is in general much smaller than 
the upper bound of Theorem 14.21 (which is probably not optimal in the case n S> d). 

A third approach, along the lines of |DS041lJ"os09] . would consist in representing tropical 
polyhedra by polyhedral complexes in the usual sense. However, an inconvenient of polyhe- 
dral complexes is that their number of vertices (called "pseudo-vertices" to avoid ambigui- 
ties) is exponential in the number of extreme rays [DS04] , Hence, the representations used 
here are more concise. This is illustrated in Figure [8] (generated using Polymake), which 
shows an intersection of 10 signed tropical hyperplanes, corresponding to the "natural" 
pattern studied in ACK0!) . There are only 24 extreme rays, but 1215 pseudo- vertices. 

Experiments. Allamigeon has implemented Algorithm ComputeExtreme in OCaml, as 
part of the "Tropical polyhedral library" (TPLib), http://penjili.org/tplib.html, Tabled] 
reports some experiments for different classes of tropical cones: samples formed by several 
cones chosen randomly (referred to as mdx where x is the size of the sample), and signed 
cyclic cones which are known to have a very large number of extreme elements [AGK09 . 
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The successive columns respectively report the dimension d, the number of constraints n, 
the size of the final set of extreme rays, the mean size of the intermediary sets, and the 
execution time T (for samples of "random" cones, we give average results). 

The result provided by ComputeExtreme does not depend on the order of the in- 
equalities in the initial system. This order may impact the size of the intermediary sets and 
subsequently the execution time. In our experiments, inequalities are dynamically ordered 
during the execution: at each step of the induction, the inequality ax < bx is chosen so 
as to minimize the number of combinations (ag^g 1 © (bg l )g J . This strategy reports better 
results than without ordering. 

We compare our algorithm with a variant using the alternative extremality criterion 
which is discussed in Sect. H] and used in the other existing implementations |CG MQ, 
IAGG08) . Its execution time T 1 is given in the seventh column. The ratio T/T' shows 
that our algorithm brings a huge breakthrough in terms of execution time. When the num- 
ber of extreme rays is of order of 10 4 , the second algorithm needs several days to terminate. 
Therefore, the comparison could not be made in practice for some cases. 

Table Q] also reports some benchmarks from applications to static analysis. The ex- 
periments oddeveni correspond to the static analysis of the odd-even sorting algorithm of i 
elements. It is a sort of worst case for our analysis. The number of variables and lines in each 
program is given in the first columns. The analyzer automatically shows that the sorting 
program returns an array in which the last (resp. first) element is the maximum (minimum) 
of the array given as input. It clearly benefits from the improvements of ComputeEx- 
treme, as shown by the ratio with the execution time X" of the previous implementation 
of the static analyzer [AGG08J. 
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